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Abstract. We analyze the role of orbital degeneracy in possible magnetic and orbital 
instabilities by solving exactly a two-site molecule with two orbitals of either e g or ti g 
symmetry at quarter-filling. As a generic feature of both models one finds that the 
spin and orbital correlations have opposite signs in the low temperature regime when 
the orbitals are degenerate, in agreement with the Goodenough-Kanamori rules. While 
Hund's exchange coupling Jh induces ferromagnetic spin correlations in both models, 
it is more efficient for t2 g orbitals where the orbital quantum number is conserved 
along the hopping processes. We show that the ground state and finite temperature 
properties may change even qualitatively with increasing Coulomb interaction when 
the crystal field splitting of the two orbitals is finite, and the Goodenough-Kanamori 
rules may not be followed. 
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1. Introduction 

The Hubbard model has been employed for a long time as the conceptually simplest 
model which might explain metallic ferromagnetism of itinerant electrons |2| and 
localization due to Coulomb repulsion [2]. Unlike initially expected, this model does 
not easily yield ferromagnetism on the hypercubic lattice and some additional conditions 
have to be satisfied to stabilize a ferromagnetic (FM) phase. As one of very few exact 
result in this many-body problem, Nagaoka established long ago [I] that in the limit of 
infinite local Coulomb interaction U = oo, a single hole doped into a half-filled system 
leads to the FM ground state. There are several indications that this remains valid for a 
finite concentration of holes jSllHI- Nevertheless the ground state of the Hubbard model 
on the square lattice may only be FM for values of U far larger than the bandwidth 

Ferromagnetism may be also promoted by a particular lattice or band structure 
jH], or even by disorder jH]. Lieb first showed that a half- filled flat band induces a 
net magnetization [Oj. Furthermore, Hirsch and others have focused on the effect of 
additional off-diagonal matrix elements of the Coulomb interaction |lUj . 

A major step towards understanding the physics of metallic ferromagnetism in 
transition metals such as Fe, Co, and Ni, was the suggestion that orbital degeneracy 
might play a crucial role. It was first pointed out by Slater, Statz, and Koster [TT] 
and then stressed by van Vleck [T2] that in the presence of degenerate orbitals, Hund's 
exchange coupling J H favours local triplet spin configurations of two electrons occupying 
different orbitals. This has important consequences for ferromagnetism as local moments 
are formed and thus Jh helps to stabilize magnetic phases, including the FM one [T3] . 
Roth examined the doubly degenerate model at quarter-filling in three dimensions [14J. 
She found that the ground state is a spin triplet and orbital singlet for two sites, i.e., the 
system forms an orbital superlattice structure in which two sublattices have different 
orbitals occupied by electrons at each of them. This interrelation between staggered 
orbital order and spin ferromagnetism was next emphasized by Kugel and Khomskii 
[T5] who derived an effective strong coupling Hamiltonian with coupled spin and orbital 
degrees of freedom, extended further by Cyrot and Lyon-Caen [TH] by on-site pair 
hopping, and by Inagaki [T7]. These seminal papers started a new field — spin-orbital 
physics in correlated transition metal oxides ^H], where superexchange models derived 
in the strong coupling regime provide a theoretical background for understanding both 
magnetic and optical properties ^H]- Numerous spin-orbital models have been derived 
in the regime of large Coulomb interactions, both with e g [20] and with t<i g EI] orbital 
degrees of freedom and are currently investigated. 

While systems of higher dimensionality are clearly the ones of most interest |2*2*| ESj , 
significant insight into the complementary behaviour of spin and orbital degrees of 
freedom was obtained in a one-dimensional (ID) model. Indeed, many of essential 
features of such ID systems were established by quantum Monte Carlo simulations [24J, 
Exact Diagonalization (ED) studies [2^1 [2^1 |^Z| , the combination of these two approaches 
|28j . and Density Matrix Renormalization Group (DMRG) method [2H]. Spin, charge, 
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and orbital correlations in the ID Hubbard model with t 2g orbitals at several densities 
have been also examined in a very recent study by Xavier et al. [HO], which combines 
DMRG and the Lanczos technique. 

Unfortunately, in a vast majority of studies, a conceptually simplified model, i.e., 
the degenerate Hubbard model with equivalent orbitals [22j or at best with different 
bandwidth obtained with diagonal hopping [23] , has been investigated so far. However, 
it has been recently shown within dynamical mean field theory that a finite on-site 
hybridization between orbitals enhances the charge and orbital fluctuations and plays 
a significant role especially for the orbital-selective Mott transition [3T]. In fact, a 
proper description of a system with e g orbitals involves a more complex kinetic energy 
Hamiltonian as in this case the orbital flavour during the intersite hopping is not 
conserved, and this is likely to lead to partial orbital polarization which is expected 
to strongly modify the magnetic instabilities |S21ISS|, resulting in a rich phase diagram. 
In addition only little is known in analytic form about these models, besides expansions 
in t around the atomic limit which result into the celebrated Goodenough-Kanamori 
rules |3U E2] and spin-orbital models. In particular cases, such as half-filling in two- 
band models, they corroborate the results obtained from weak coupling approaches, 
and therefore gaining further qualitative insight into the corresponding problems is 
unlikely. In contrast, at quarter filling, weak coupling approaches are at odd with the 
strong coupling expansions [22!, an d the situation remains controversial. Fortunately 
this situation can be studied rigorously on the analytical level, both for e g and t 2g 
orbitals, in a two-site molecule. 

The paper is organized as follows. Firstly, we introduce a realistic model with 
two e g orbitals and compared it with the one for t 2g orbitals in section [21 Secondly, 
the solutions for both models are given and compared with each other in section El 
Thereby, we identify characteristic differences in the behaviour of both orbital degrees 
of freedom [SB!- Furthermore, in section El we verify whether the phenomenon of 
a complementary behaviour of the spin and orbital flavours is also a characteristic 
feature of this particular case by analyzing spin and pseudospin correlation functions. 
Comparisons to strong coupling expansions are presented as well. Next, in section 0] we 
investigate the influence of finite crystal field splitting on the ground state properties. 
Finally, section El summarizes the paper and gives general conclusions. 

2. Model Hamiltonian for e g and for t 2g orbitals 

The magnetic and orbital instabilities within the e g band become especially relevant in 
the context of doped nickelates La2- a Sr x Ni04 (LSNO), where interesting novel phases 
including the stripe order were discovered [21]. Even though doped nickelate LSNO 
is isostructural with its cuprate counterpart La 2 - a; Sr :r Cu04, its electronic degrees of 
freedom are more complicated. In fact, a realistic Hamiltonian for LSNO must contain, 
besides the x 2 — y 2 orbital states deciding about the properties of the cuprates, also the 
3z 2 — r 2 orbital states, so as to account for the actual filling with two holes and for the 
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high-spin state (S = 1) in the stoichiometric compound. Such a model of interacting e g 
electrons in (a, b) plane may be written as follows, 

H = H kin + Hi nt + H c f, (1) 

with two orbital flavours: \x) ~ x 2 — y 2 and \z) ~ 3z 2 — r 2 forming a basis in the orbital 
space. The kinetic energy is described by, 

Hkin = ^^2t i fc\ aa C jf3a , tff = -- ( 1 )' ^ 

{ij) af3a \ / 

where t stands for an effective (dda) hopping matrix element due to the hybridization 
with oxygen orbitals on Ni— O— Ni bonds, and the off-diagonal hopping t x ^ along a and 
b axis depends on the phase of the \x) orbital along the considered cubic direction. 

For our purpose it is most convenient to consider an e g orbital basis consisting of 
a directional orbital |£) along the molecular bond and a planar orbital |£) orthogonal 
to \() |28j . Pairs of such orthogonal orbitals defining a new basis might be obtained by 
the following transformation of the original orbital basis {\x), \ z)}, 



10 \ I cos i sin f \ ( l*> 

10 ) \ -sinf cosf ) \ \x) 



(3) 



with the angle 9 = ±2n/3 depending on whether one considers the bond along a or b 
axis. This rotation leads to a simple diagonal form of the hopping matrix, 



allowing only for intersite transitions between two directional \() orbitals along each 
bond We compare this case with a frequently studied degenerate Hubbard model 
with equivalent orbitals. Such a model describes the dynamics of two active t2 g orbitals, 
e.g., for a molecular bond along a axis, \() ~ \xz) and |£) ~ \xy), so that the hopping 
matrix is diagonal, 

The electron-electron interactions are described by the on-site terms, which we 
write in the following form j39j . 

H in t = U ^2(n ix1 n ixl + n iz1 n izl ) + (U - \J H ) ^ n ix n iz 

i i 

— 2Jh 2j Sj x • S iz + Jh S ^{ c \x'\ C \xl C izi C iz] + Az\Az{ C ixi C ix\) ■ (6) 
i i 

Here U and Jh stand for the intraorbital Coulomb and Hund's exchange elements, 
whereas n ia = n iaa is the electron density at site i in a = x, z orbital state. 

The last term H c f stands for the uniform crystal- field splitting between \x) and \z) 
orbitals along the c axis, 

Hcf — \Eq ^2( Uixa ~ n iz°)- (^) 
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The splitting of two e g orbitals occurs due to the tetragonal Jahn- Teller distortion 
of the NiOg octahedron. In La2Ni04, however, the octahedron, with the Ni— O— Ni 
bond lengths being 1.95 (2.26) A in-plane (out-of-plane) |10], respectively, is much less 
distorted than the CuOg octahedron with 1.89 and 2.43 A bond lengths in La2Cu04 
[4~T] . which reflects the difference in electron filling. In what follows we consider only a 
realistic positive E favouring, due to elongated octahedra, the \z) occupancy over the 
\x) occupancy by the e g electrons in doped compounds. 



3. Physical properties at orbital degeneracy 

3.1. Classification of eigenstates using the symmetry properties 

In this section we present an exact solution of a two-site molecule at quarter filling with 
two degenerate orbitals of either e g or t2 9 symmetry. Although it is straightforward 
to solve the present problem numerically, it is more instructive to find the solution 
analytically. Along this process several important aspects will be clarified. 

While the total spin operator commutes with the Hamiltonian and its z-component 
can be used to label the eigenstates this does not hold for the total orbital pseudospin 
operator T = ^ Tj, where its three components are given by, 

a a a 

In order to distinguish the third component of the pseudospin from that conventionally 
used for a bond along c axis (% z ), we have labelled it with the ( index. Note however 
that in contrast to the spin operator, the ^-component of the pseudospin operator, 

^ C = E^' ( 9 ) 

i 

does not commute with the Hamiltonian (Q), so the states with different values of this 
observable mix with each other. Nevertheless, we will use its eigenvalues together with 
the ^-component S z of the total spin operator S to specify multiparticle states in terms 
of which we write the Hamiltonian. It is straightforward to construct explicitly all 28 
states; they are listed in |Appendix A 



In the high-spin subspaces S z — ±1, the Hamiltonian Q is decomposed into a zero 
2x2 matrix H 0a involving the two \^ aa ) states with T*> = ±1, and two 2x2 matrices: 

^ = (<*a<*,vi) i , \ ;r ;r ) i i m 




and, 



Now we are left with the S z = subspace. In the S = 1 sector one recovers the same 
eigenenergies: two zeros corresponding to the above = ±1 localized states with 
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= ±1, and the eigenvalues following from the matrices Hi a , and H 2a - In the S = 
subspace, using the states with = ±1, the following Hamiltonian matrices are found 
in addition: 



u 



( 2t 

2t cc 


V 



cc 





Jh 



\ 





J H U 2t 



2t 



/ 



m) 
v i*r> J 



(12) 



and, 



:^ 2 u <^3 




Finally, in the sector with = 0, the Hamiltonian matrices read, 

l*s) 




U — Jh ttf + % 
t cc + % 



Eigenenergies for t2 9 orbitals 



Jh t 



CC 



l*F> 



(13) 

(14) 
(15) 



Let us consider first a model with two active and equivalent ti g orbitals. The third 
orbital may be neglected in various contexts: {%) either in a d 1 configuration when a 
crystal field raises the energy of the third orbital above the other two we consider, or (ii) 
in a d 3 low spin configuration when a crystal field lowers the energy of the third orbital 
below the other two (the former one being filled), or (Hi) in a d 5 low spin configuration 
when a crystal field lowers the energy of the third orbital below the other two in which 
case one meets the problem of one hole in two orbitals. From the few layered perovskites 
that may correspond to these cases let us mention NdSrCr04 [12] ■ We note that the 
distortions in layered perovskites differ from the ones in cubic perokskites, where they 
are typically trigonal. For this model every matrix in equations ()12)) - fT5j) can be easily 
diagonalized and the corresponding eigenvalues are listed in the | Appendix B For clarity 
we have used the symmetry and classified them into triplet (S = 1) and singlet (S = 0) 
subspaces. 

To get more insight into the competition between the tendencies towards the 
antiferromagnetic (AF) and FM ground states, we now discuss the lowest energy 
eigenstates in the strong coupling (large U) limit. As expected, the lowest energy, 

E T 3 =i(U-3J B - V(tf-3J H ) 2 + 16t 2 ) , (16) 

obtained by the diagonalization of the matrix in equation ([lip, belongs to spin triplet 
coexisting with the pseudospin singlet. In the case of strong on-site interorbital repulsion 
(U — 3Jh t), the lowest high-spin energy (fTB^) reads, 

At 2 
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Table 1. Lowest eigenenergies of the model for ti g orbitals at U = 8t, obtained for 
the representative values of Jjj = 0, U/8 and E//4. The eigenstates \n) are specified by 
the total spin S n and the expectation value of the ^-component of the total pseudospin 
JHJ). Triplet states with S n = 1 have three components S* = ±1,0. 



Jh = 




Jh = 


U/8 




Jh = 


U/4 




E n /t S n 


n 


E n /t 


S n J 


rC 

' n 


Ejt 


S n 


n 


0.4721 1 





-0.7016 


1 





-1.2361 


1 








±1 


-0.5311 





Q 


-0.6056 































Note however that finite Jh could significantly reduce the value of the interorbital 
repulsion U — 3 Jh, so that it would no longer be much larger than t. As a consequence, 
significant corrections to the above result obtained to second order in t/(U — 2>Jh) are 
expected in this case. Analogously, the lowermost low-spin energy results from the 
matrix in equation ()12j) . which is degenerate with the one obtained from the matrix in 
equation (JHJ), 

/ \ At 2 

41 9 = -\(u-j h - y/(u-j H )* + m*) ~ -7737^ ( 18 ) 

in the strong coupling regime. 

Comparison of equation (j!8|) with the lowest high-spin energy given by equation (|17|) 
allows one to draw interesting conclusions about conditions required for ferromagnetism. 
It is apparent that and E^ are degenerate for Jh = 0. However, even 

infinitesimally small Jh > lifts this degeneracy and might give rise to spin 
ferromagnetism combined with the alternating orbital (AO) pseudospin correlation. 

As an illustrative example, we present in table Q the lowest eigenenergies E n 
obtained from the ED of the Hamiltonian ([!} for the t2 9 orbitals in the strong coupling 
regime U = 8t for a few values of Jh- We have specified them in terms of the total spin 
S n = 0,1. Another quantity used for the classification of states is the (^-component of 
the total pseudospin operator, see equation Q. However, even in this case it is not a 
good quantum number being modified by the 'pair-hopping' processes from one orbital 
to the other present in equation (jHJ. Indeed, in matrix ([12)1. a sector consisting of |<&]~) 
and I $2) states, carrying the = —1 pseudospin flavour, is coupled to the one in terms 
of I $3 ) and I $4 ) states, carrying the opposite T"> = 1 flavour. Similarly, sectors with 
different pseudospin flavours are mixed in matrix (fTo^ . Thus, one has to determine the 
expectation values of by a direct evaluation using the eigenstates of the Hamiltonian. 

The ground state of the t2 9 model with finite Jh is a spin triplet accompanied by 
a pseudospin singlet (cf. tableHJ. Notice however that its energy — 1.2361i obtained 
for Jh = U/4 differs vastly from the value E^ = —2t estimated roughly from equation 
(fTTj) . The reason of this discrepancy is the failure of the second order perturbation 
theory controlled by t 2 /{U — 3Jh), being here of order 0(t). Therefore, it can be used 
only for qualitative arguments in this regime, while it works better for smaller Jh = U/8, 
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Table 2. The same as in tabled but for e g orbitals, and four representative values of 
Hund's exchange: J H = 0, 17/8, U/A, and U/3. 



Jh = 




Jh = 


[7/8 


Jh 


= U/4 


Jh 


= U/3 




E n /t S n 


n 


E n /t S n 


T C 

n 


E n /t 


S 


E n jt 


S n 


■* n 


0.4721 


-l 


-0.4784 


-0.9985 


-0.4983 


-0.9935 


-1.0 


1 


0.0 


0.1231 1 





-0.1926 1 


0.0 


-0.4142 


1 0.0 




1 


0.0 


1 





1 


0.0 




1 0.0 


-0.5204 


- 


0.9877 








-0.1401 


0.0 


-0.1623 


0.0 


-0.1813 





0.0 











0.0 




0.0 







0.0 



yielding there — 0.8i, a value much closer to the exact energy — 0.7016t. 
3.3. Eigenenergies for e g orbitals 

Turning now to the e g orbital model and recollecting the immobile |£) orbital flavour 
with tg£ = 0, one immediately notices two sets of two identical 2x2 subspaces spanned 
by the = states, described by the matrices in equations (|TIj) - ([T5|l . which readily 
yield two doubly degenerate singlets (c/. | Appendix Q . Diagonalizing the submatrix (fT^j) 
yields the three eigenvalues given in |Appendix C| The first of them (jC.3|) corresponds 
to the lowest energy low-spin state. In the strong coupling regime, i.e, for (U — Jh) t, 
one finds 

At 2 U , x 

K - -JJT^- («•) 

It has to be compared with the lowest eigenvalue of the S z = 1 subspace, the lowest 
energy doubly degenerate spin triplet, 

E t = l(u-3J H - ^(U -3J H ) 2 + 4t^, (20) 

obtained by the diagonalization of the matrices (|Tnj) and pi)) . It corresponds, as in 
the model with t 2g orbitals, to the pseudospin singlet (T n = 0). In the strong coupling 
regime it yields the lowest high-spin energy, 

E e T 9 ~ . (21) 

Note that the energy is much higher than that for the t 2g model, see equation (|T7)l . as 
only one electron is mobile. 

Based on the lowest energy excitations (|19|) and (|21j) one can easily notice a striking 
difference between t 2g and e g orbitals with respect to the ground state. Indeed, in 
contrast to the t 2g model where even infinitesimally small Jh > lifts the degeneracy 
between the lowest energy singlet and triplet excitations and stabilizes the FM spin 
correlation, one can expect that the singlet state with the energy Eg 3 ()C3|) remains the 
e g ground state up to Jh ~ U / 4. 

The corresponding exact eigenenergies of the Hamiltonian (P) with e g orbitals 
obtained in the strong coupling regime U = 8t for a few values of Jh are listed in 
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table 121 The data show that increasing Jh diminishes the energy difference between the 
lowest singlet and the first two triplet excited states. Remarkably, however, even for 
unrealistically large Jh = U/4, the e g ground state is still a singlet with almost fully 
occupied mobile pseudospin \() orbitals (cf. table |2J) contradicting our predictions from 
the strong coupling regime. It follows from the approximate low-spin state energy E s " 
(|T9*j) which overestimates the tendency towards ferromagnetism due to the performed 
expansion, and the transition from the singlet to the high-spin state occurs only for 
J H ~ 0.27C7. Indeed, the ground state found for Jh = U/3 is a triplet, with the energy 
well below that of the singlet state (see table EJ). Therefore, although Hund's exchange 
coupling Jh is a driving force of the FM spin correlations in both models, it is decisively 
more efficient in the case of t 2g orbitals, where infinitesimal Jh stabilizes the high-spin 
state as both pseudospins are mobile. 

Actually, one can rewrite the interaction term of the Hamiltonian Hi nt © as a 
superposition of interactions in different channels, involving the density, magnetization, 
orbital polarization, magnetic orbital polarization, and on-site orbital flip term, 
respectively [32] ■ I n each of them the interaction strength is different: it is systematically 
attractive (repulsive) in the magnetic (density) channels, while it turns repulsive in the 
orbital polarization channel for 5J# > U. As a result, one expects a transition in the 
ground state from pseudospin triplet to singlet with increasing Jh/U . This is precisely 
the transition we discussed above for e g electrons. In contrast, the magnetic instability 
takes over for t 2g electrons, and the spin triplet ground state is found at any J H > 0. 

3.4- Correlation functions and susceptibility at orbital degeneracy 

We now turn to the temperature dependence of the on-site (S 2 ) and intersite (Si • S 2 ) 
spin correlations, and the on-site (T 2 ) and intersite (Ti • T 2 ) pseudospin correlations. 
The latter yield information about an orbital state together with orbital correlation 
between neighbouring sites. 

In figure [T] we present the temperature dependence of the spin and pseudospin 
correlation functions as well as both susceptibilities of the model with e g orbitals. We 
have set Hund's exchange coupling to be Jh/U =1/8 (dashed line) and Jh/U = 1/4 
(solid line). At low temperature, one expects charge localization as the system is in the 
strong coupling regime U = St. Indeed, the local spin moment (S 2 ) reaches virtually 
the magnitude S(S + 1) = 3/4 for the spin S = 1/2. A rise of (S 2 ) above this value 
upon increasing temperature is caused by thermal excitations to triplet states. They are 
favored by Hund's interaction and form local high-spin configurations. Consequently, 
the rise of (S 2 ) is larger for stronger Jh = U/A. Next, the intersite spin correlation 
function (Si ■ S 2 ) indicates the low-spin (AF) nature of the ground state, whereas the 
corresponding pseudospin function (Ti-T 2 ) illustrates the ferro orbital (FO) pseudospin 
correlation. 

The results obtained within the t 2g model are qualitatively different. First, at 
Jh = 0, apart from spin SU(2) symmetry there is an additional SU(2) symmetry for 
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Figure 1. Temperature dependence of local moments (a) spin (S^), (b) pseudospin 
{T i ), and of intersite correlations (c) spin (Si -82), (d) pseudospin (TVT 2 ), as obtained 
for the model with e g orbitals, for increasing values of Hund's exchange: Jh = U/8 
(dashed line), Jh = U/4 (solid line) and Jjj = [7/3 (dotted line). Parameters: U = 8t, 
-Bo -0. 



orbital degrees of freedom, resulting in higher SU(4) symmetry [131 an d intersite spin 
and pseudospin correlation functions are both negative and identical. Second, at finite 
Jh positive (Si • S 2 ) indicates the FM nature of the ground state supported by the 
pseudospin singlet with negative (Ti • T 2 ) correlations shown in figure El 

The gradual increase of triplet intersite correlations (Si • S 2 ) observed for the e g 
model is also well recognized in the spin susceptibility, see figure E^a). Upon taking the 
logarithm of \s we find a typical AF behaviour with a characteristic cusp at the crossover 
temperature T c (e g ). Obviously, the AF array that sets in has zero net magnetic moment 
at the temperature below T c (e 9 ) and this explains the observed cusp in Xs- Only for 
large Jh = U/3 the ground state is a triplet, and the spin susceptibility changes to 
the Curie- Weiss type. In contrast, in case of t 2g model [figure Of b)], Xs is large at low 
temperature already for J H = as the ground state has degenerate singlet and triplet 
states. Increasing Jh gives here only moderate enhancement of the spin susceptibility 
as the degeneracy of the ground state is then removed. 



3.5. Specific heat and entropy at orbital degeneracy 

Different energy spectra of the e g and t 2g systems result in quite different temperature 
behaviour of the specific heat C and the entropy S, as shown in figures |U and Consider 
first the e g system. As depicted in figure IH a low temperature peak of the specific heat 
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Figure 2. The same as in figure^but for the model with t^g orbitals, and for: Jh = 
(dotted line), Jh = U/8 (dashed line), Jh = U/A (solid line). 



coincides with the characteristic kink in the susceptibility Xs [c/. figureHHa)]. Comparing 
the position of the peaks corresponding to Jh = U/8 and Jh = U/4, one finds that 
increasing Jh reduces T c (e g ), as the crossover to the high-spin state is approached. Note, 
however, that for Jh = U/A the low temperature peak in the specific heat splits into two. 
The first one corresponds to a transition from the singlet ground state to the first two 
triplet excited states with the excitation energy AEi/t = 0.0841, whereas the second 




Figure 3. Temperature dependence of the spin susceptibility log(xs/7 2 A0 per site 
(N = 2) as obtained for two values of Hund's exchange, Jh — U/8 (dashed lines) and 
Jh = U/A (solid lines), for: (a) e g model, and (b) tig model. For comparison the 
susceptibilities obtained with Jh = U/3 for the e g model and with Jh = for the ti g 
model are shown by dotted lines. Parameters: U = 8t, Eq = 0. 
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Figure 4. Temperature dependence of (a) the specific heat C per site and (b) the 
entropy S, as obtained for e g orbitals with: Jh = U/8 (dashed line), Jh = U/A (solid 
line), and Jh — U/3 (dotted line). Parameters: U — 8t, Eq = 0. 



one appears due to the higher energy charge excitations, with the excitation energies 
AE 2 /t = 0.336 and AE 3 /t = 0.498, respectively. In contrast, when Jh = U/8 the 
excitation energy into the first excited state is much larger AE[/t = 0.286, whereas the 
other excitation energies are nearly unaltered: AE' 2 /t = 0.338 and AE' 3 /t = 0.478. This 
results in a single broad low temperature peak. High temperature peaks occur due to 
thermal excitations which create double occupancies and lead to charge derealization, 
well seen in the suppression of (S?). Finally, at large Jh = U/3 all peaks in the specific 
heat merge into a single broad peak at higher temperature. 

Turning now to the temperature dependence of the entropy depicted in figure Efb), 
one finds S ~ in the low T regime for the singlet ground state. Basically, the overall 
rapid increase of the entropy around T c (e g ) in figure Efb) is very much the same for 
Jh = U/8 and Jh = Z7/4, corresponding to the low temperature peak in the specific 
heat. However, a detailed behaviour of S depends on Jh- In the regime of small 
Hund's exchange Jh < U/8, where Ls = 16 singly occupied states are well separated 
from doubly occupied ones, S possesses a point of inflection S = k B In 16 at k B T ~ t, 
which follows from the gapped character of the specific heat. In contrast, for the larger 
Jh = U/4, the entropy increase starts at lower temperature as the energy of the singlet- 
triplet excitation is low. Here the gap between singly and doubly occupied states is 
smaller and the corresponding point of inflection is less transparent. The limiting value 
S = /c£m28 results from the calculation performed in the canonical ensemble. 

The situation is quite different in the case of the t2 9 model. Unlike the e g case, 
increasing J H shifts T c (t 2g ) towards higher temperatures (figure EJ). As a result, the 
crossover temperatures of both systems differ substantially, especially in the large 
Jh = U/A regime. Indeed, from the position of the low temperature peak of the 
specific heat in figures Ufa) and|Sfa), one can read off that k B T c (e g ) = 0.025t, whereas 
kBT c (t2 9 ) = 0.35t. The origin of this marked difference is certainly the fact that 
increasing J H diminishes (enlarges) the gap between the spin singlet (triplet) ground 
state and the first excited triplet (singlet) state of the e g fag) system, respectively (cf. 
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Figure 6. Zero temperature intersite correlation functions: spin (Si • S2) (solid lines) 
and pseudospin (Ti • T 2 ) (dashed lines) as functions of the Stoner parameter U + Jh 
in the two-band model at Jh — U/A for: (a) e g orbitals, and (b) ti g orbitals. 



tables H] and EJ). In contrast to the e g case, the entropy is finite in the low temperature 
regime, and approaches either the value /celn3 for both Jh = U/8 and Jh = U/A or 
k B \n.Q when J H = 0, see figure E(b). Indeed, for J H = 0, the ground state corresponds 
to six degenerate states — three of them constitute a spin triplet, whereas the others 
are singlets, as shown in table [TJ However, any finite Jh > splits up these states and 
leads to the triplet ground state with the entropy S = ks In 3. As in the e g case, S has 
a clear point of inflection only in the small Jh < U/8 regime. 

To summarize this section, we compare in figure |H1 the intersite spin (Si • S2) and 
pseudospin (Tj ■ T 2 ) correlation functions for the e g system [cf. figure EJa)] an d for the 
t 2g system [cf. figure Efb)] as a functions of the Stoner parameter U + Jh- Note that 
(Si • S2) is finite and negative even in the noninteracting U = limit due to the Pauli 
principle. On the one hand, the e g results illustrate the AF correlations in the ground 
state (negative (Si ■ S 2 )) owing to the preferred mobile pseudospin \Q. On the other 
hand, the ground state of the ti g system with a finite interaction is a spin triplet as 
(Si • S2) is positive and a pseudospin singlet as the corresponding correlation function 
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Table 3. Eigenenergies of model with t^g orbitals as obtained at U — 8t with 
a finite crystal field Eq — It acting perpendicular to the chain, for two values of 
Hund's exchange: Jh = U/8 and Jh = U/4. The spin quantum number S n and the 
expectation value of T% operator JjJJl in each eigenstate \n) are given. 
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(T 1 ■ T 2 ) is negative. We therefore conclude that the ground state properties strongly 
depend on the orbital correlations, in agreement with the Goodenough-Kanamori rules 



4. New features at finite crystal field 

4-1- Correlation functions and susceptibility at finite crystal field 

We now discuss the influence of the crystal field splitting given by equation ((7|). Here we 
are interested in the nontrivial case of crystal field acting along the c axis perpendicular 
to the chain. Hence, as we have been working with the e g orbital basis consisting of a 
directional orbital along the molecular bond \() and an orthogonal to it planar orbital 
one needs to rotate the field (JJJ) expressed in the original orbital basis {\x), \z)} by 
the same angle 9 = 2tt/3 which enabled us to simplify the form of the hopping matrix 
(J2J) into P|). Making the inverse transformation to (J3J) in equation Q, one obtains the 
crystal field term, 

H <* = k E o Yl [ cos d ( c h* c ita ~ 4<;* c i(v) + sin 8( c lt* c ic* + 4i* c iZ*) ■ ( 22 ) 

ia 

Note that cos# is negative so the field fl2H) favours the |^) occupancy over the \() 
occupancy, as it should. 
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Table 4. The same as in table OH but for e g orbitals. 
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Unfortunately, except for the Hamiltonian matrix (Jll)) . the form of the other 
matrices is now considerably more involved due to offdiagonal elements in the crystal 
field in equation (f2~2"|) which couple states with different T**. In general, it is not possible 
to obtain analytic expressions for the eigenvalues and one has to resort to a numerical 
diagonalization. However, due to the equivalent hopping amplitudes, the eigenvalues 
of the t 2g model should be independent of the rotation angle 9, i.e., one has to get the 
same energy spectrum for the field (}2~2*|) with finite 9, as well as for the diagonal in 
pseudospin space field of the form E T^. It acts along the chain and corresponds to 
9 = in equation (|22jl. 

Numerical values of the eigenvalues obtained for t2 g orbitals for Eq = It and 
9 = 2n/3 are listed in table 01 Although the energy spectrum of the t 2g model indeed 
does not depend on the field direction, expectation values of in the Hamiltonian 
eigenstates certainly do, as a finite rotation angle 9 enables the mixing of states with 
different values of T^. For example, the initial pseudospin = ±1 of spin triplets 
with the energy Eo/t = ±2 is reduced up to = ±0.5 by the field with 9 = 27r/3 (cf. 
table El), whereas it is conserved when 9 = 0. 

Contrary to the -E^ = case with a spin triplet and pseudospin singlet as the 
t 2g ground state at finite Jh > 0, finite positive Eq suppresses the AO pseudospin 
correlation and stabilizes a spin singlet with a positive value of ~ 0.5, inducing FO 
correlations. Remarkably, the effect of the crystal field on the e g ground state is just 
the opposite, as reported in table EJ Namely, by lifting the degeneracy of pseudospin 
flavours it promotes the immobile |£) one. Consequently, there is not that much kinetic 
energy to be gained and the Coulomb interactions start to be crucial. However, they 
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Figure 7. Temperature dependence of local moments (a) spin (S?) and (b) pseudospin 
(T?), and ofintersite correlation functions (c) spin (Si«S2) and (d) pseudospin (Ti-T 2 ), 
as obtained for the model with e g orbitals for: Jh — U/8 (dashed lines) and Jh = U /A 
(solid lines). Parameters: U = 8t and Eq = 2t. 



are noticeably better optimized by the FM spin correlation. Indeed, from table |U one 
sees a strong competition between the lowest singlet and triplet states, both with a 
positive but still smaller than in the t% g case value of ~ 0.45. Note also that it 
becomes energetically advantageous to have the spin triplet as the ground state for 
large Jh = U/A while a smaller Hund's exchange coupling Jh = U/8 drives the system 
towards the singlet in the ground state. 

The above ground states and the excitation spectra obtained at finite crystal field 
Eq = 2t influence intersite spin and orbital correlations, shown in figure [H The ground 
state of the e g system depends on the value of Jh, being a spin singlet for J H = U/8 
which results in negative (Si • S2) and a spin triplet for Jh = U/A yielding positive 
(Si • S 2 ) in the low temperature regime. In contrast, (T x ■ T 2 ) is positive (FO pseudospin 
correlations) and almost insensitive to the value of Jh at low temperature. For increasing 
temperature one finds two transitions — firstly, the intersite spin correlations weaken, 
and secondly pseudospin correlations weaken at much higher temperature, when charge 
excitations are thermally activated. This separation of the energy scales for spin and 
orbital excitations occurs both for lower J H = U/8 and higher J H = U/A. 

A finite crystal field affects drastically the behaviour of the t 2 g correlation functions 
as well (figure EJ). At low temperature negative (Si -82) reveals the AF coupling between 
spins, whereas positive (Ti • T 2 ) indicates the FO pseudospin correlation, regardless of 
Jh- In contrast to the e g case, however, spin and pseudospin correlations are reduced 
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and vanish simultaneously, indicating that spin-orbital degrees of freedom are stronger 
coupled in this case [2T] . 

The change of the magnetic correlations with increasing Jh in the e g model is 
reflected in the crossover from the AF behaviour of Xs at Jh = U/8 to the Curie- Weiss 
behaviour of Xs obtained for Jh = U/4, see figure [Ufa). In case of ti g model [figure E(b)] 
one finds the AF character of Xs in the interesting range of Jh- Altogether, as for 
E z — 0, the spin susceptibility Xs exhibits again the opposite behavior for both types 
of orbitals — the high-spin state is suppressed in the ti g model, while it can be selected 
by the crystal field in the e g model. The maximum of Xs occurs at higher temperature 
in tig than in e g model at Jh = U/8. 



4-2. Specific heat and entropy at finite crystal field 

As expected from the above results, a finite crystal field also modifies the temperature 
evolution of the specific heat C and the entropy S for e g orbitals, see figure figure E3 
The position of the high temperature peak of the specific heat at ksT ~ t is almost the 
same as the position of the strong anomalies of both on-site correlation functions (cf. 
figures [TUfa) and[7|). Depending on J H , the low temperature entropy S of the e g system 
either vanishes (for the singlet ground state as Jh < U/8) or approaches the value 
fcehiS (for the triplet ground state at Jh = U/4). Nevertheless, owing to the vanishing 
specific heat, all the curves in figure ITUT b) have a point of inflection S = /c£ln4 at 
fceT ~ O.lt. Note that in contrast to the case without a crystal field, there is no such a 
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Figure 8. The same as in figured but for the model with t^g orbitals. Dashed (solid) 
line corresponds to Jh =U/8 (Jh = U/A), respectively. 




Figure 9. Temperature dependence of the spin susceptibility \og(xsh 2 N) per site 
(N = 2) for: (a) e g model, and (b) t 2g model, as obtained for two values of Hund's 
exchange: Jh — U/8 (dashed lines) and Jh = U/4 (solid lines). Parameters: U = 8t, 
E Q = 2t. 



point when S = /cglnlG. Namely, by promoting one pseudospin over the other one, a 
finite crystal field markedly lowers the states with double occupancies, hence the usual 
gap between the singly and doubly occupied states vanishes. 

Consider now the temperature behaviour of the specific heat for t2 9 orbitals shown 
in figure HIT a). Instead of the high and low temperature peaks of the specific heat one 
observes two slightly separated peaks for J H = U/8 which merge into a wide peak at 
Jh = U/A. This confirms that spin-orbital degrees of freedom are strongly coupled as 
spin and orbital intersite correlations change at the same temperature (see figure [SJ). 
Finally, as shown in figure fTTT b). the entropy of the ti g system is almost independent of 
Jh- We note that S is entirely suppressed at low T due to a spin singlet ground state; it 
rises rapidly at k B T ~ O.lt, and approaches eventually the limiting value S = fceln28, 
not having any point of inflection. Such a behaviour is a direct consequence of a single 
broad peak in the specific heat and no separation of the energy scales for spin and orbital 
excitations. 

We close this section with a short discussion of the intersite spin (Si • S 2 ) and 
pseudospin (T x • T 2 ) correlations for the e g system depicted in figure H2fa) as functions 
of the Stoner parameter, I = U + Jh, in the presence of a finite crystal field Eq = It. 
As expected, the results illustrate the AF correlation between spins on two sites in the 
weak coupling regime U + Jh^, 7t. However, further increase of the interaction strength 
changes gradually the AF coupling into a FM one, with the immobile |£) pseudospin 
component preferred by the crystal field. Remarkably, such a transition at U + Jh — 16t 
in the presence of the same crystal field Eq = It from the AF phase into the FM one has 
been obtained on the infinite lattice in the Hartree approximation recently [^2] ■ (The MF 
phase diagram appears more involved |32j . since FM order is stabilized for 5 < I/t < 6.5, 
in contrast to what is obtained in the present study. Nevertheless this mostly reflects 
the strong competition between the FM and AF phases, since their energy do not differ 
by more than 1% in this regime [H].) Consequently, owing to a strong competition 
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Figure 10. Temperature dependence of: (a) the specific heat C, and (b) the entropy S 
for e g orbitals. Dashed (solid) lines corresponds to Jh = U/8 (Jh = U/4), respectively. 
Parameters: U — 8t and E = 2t. 




Figure 11. The same as in figure ^3 but for ti g orbitals. Dashed (solid) line 
corresponds to Jh = U/8 (Jh = U/4), respectively. 



between the singlet and triplet states with the lowest energies (c/. table EJ), respectively, 
one can conjecture that fluctuations clearly affect the intersite spin correlations. Indeed, 
comparison of (Si ■ S2) and (Si • 82)0 determined in the Ising limit (i.e., taking only 
Sf operators and neglecting the last term in the Hamiltonian ©) shows that dynamics 
acts to reduce the AF coupling between spins. In contrast, finite crystal field affects 
(Ti ■ T 2 ) only slightly, and the pseudospin correlations are always positive for increasing 
/ regardless of E , which implies the FO coupling between two pseudospins [cf. figures 
Eta) and 112( a)]. Thus, in the regime of large Stoner parameter / one finds that the 
classical Goodenough-Kanamori rules jSl] are violated as the crystal field stabilizes a 
particular orbital configuration. In contrast to quantum spin-orbital entanglement |35j . 
this effect here is static and may play an important role for the observed magnetic and 
orbital correlations in transition metal oxides, as discussed recently on the example of 
LiNiOa EE]. 

Turning now to the t2 g model with positive (Si • S2) and negative (Ti • T 2 ) in the 
degenerate case shown in figure Efb), the situation is also changed drastically by a finite 




Figure 12. Zero temperature intersite correlation functions: spin (Si -82) (solid lines) 
and pseudospin (Ti • T 2 ) (dashed lines) as functions of the Stoner parameter U + Jh 
for: (a) e g orbitals, and (b) t 2g orbitals. Parameters: E — 2t, Jh — U/A. Dotted 
line shows the spin-spin correlation function (Si • 82)0 obtained in the Ising limit (see 
text). 



crystal field as depicted in figure IT2T b). Indeed, the resulting (Si • S2) is then negative 
revealing the AF nature of the ground state. However, due to the energy gap between a 
singlet ground state and the first excited triplet state (cf. table EJ), fluctuations modify 
the value of the intersite spin correlations only slightly, so that (Si • S2) and (Si • 82)0 
almost overlap. Finally, positive (Ti ■ T 2 ) shows that the AO pseudospin correlations 
found before at E = are suppressed. 

5. Summary 

We have obtained the exact numerical results for a two-site Hubbard model with two e g 
orbitals at quarter-filling. By an appropriate transformation of the original orbital basis 
{\x), \z)}, into a basis consisting of a directional orbital \() along the molecular bond 
and a planar and orthogonal to it orbital we have simplified the hopping matrix 
making use of the properties of e g orbitals — the electrons in \() directional orbitals are 
mobile, while those in orthogonal |£) orbitals are fully localized [SB]- The results were 
compared with the doubly degenerate Hubbard model with two equivalent t 2g orbitals 
active along the molecular bond direction. 

As the first important result, a striking difference between t 2g and e g orbitals with 
respect to the ground state has been established. Indeed, in contrast to the t 2g model, 
where even infinitesimally small Hund's exchange Jh > lifts the degeneracy of the 
lowest energy singlet and triplet states and stabilizes the FM spin correlation, we have 
found that a spin singlet ground state survives up to Jh — 0.27U for e g electrons. In this 
regime of parameters, the intersite spin correlations (Si • S2) indicate the AF nature of 
the ground state, whereas the pseudospin function (Ti-T 2 ) illustrates the FO pseudospin 
correlation. On the contrary, in the t 2g model, except for Jh = when the intersite 
spin and pseudospin correlation functions overlap and are negative, positive (Si • S2) 
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demonstrates the FM nature of the ground state supported by the pseudospin singlet 
with negative (T x ■ T 2 ), i.e. AO correlations. Therefore, a complementary behaviour of 
the spin and orbital flavours is observed as a generic feature of both models, in agreement 
with the Goodenough-Kanamori rules 

We have further demonstrated the influence of a finite crystal field Eq on the ground 
state. On the one hand, it suppresses the FM spin correlation in the t 2g model and 
stabilizes the spin singlet ground state with negative (Si ■ S 2 ), accompanied by positive 
(Ti • T 2 ). On the other hand, the effect of the crystal field on the e g ground state is the 
opposite one. Namely, by lifting the degeneracy of the pseudospin flavours it promotes 
the immobile |£) one. Consequently, not much kinetic energy can be gained and the 
Coulomb interactions start to dominate. However, they are noticeably better optimized 
by the FM spin correlation. Therefore, it becomes energetically advantageous to have 
the spin triplet and the ground state for large Jh = U/4 yields positive (Si • S 2 ), while 
a smaller Hund's exchange coupling J H = U/8 drives the system towards the singlet in 
the ground state which results in negative (Si • S 2 ). In contrast, (Ti ■ T 2 ) is positive 
and almost insensitive to the actual value of Jh in the range < Jh < U/4. 

Next, we have investigated how the e g ground state properties evolve as a function 
of the Stoner parameter I = U + Jh in the presence of a finite crystal field E = 2t. 
As expected, we have found the AF correlation between spins on two sites in the weak 
coupling regime U + Jh ^ 7t. However, further increase of the interaction strength 
changes gradually the AF coupling into a FM one, with an immobile |£) pseudospin, 
preferred by the crystal field. Finally, by comparing the correlation function (Si-S 2 ) with 
the one determined in the Ising limit (Si-S 2 ) , we have elucidated the role of fluctuations 
in the intersite spin correlation function and have observed that the dynamics helps to 
reduce strongly the AF coupling between spins in this case. In contrast, when the ground 
state is FM, fluctuations only slightly modify the value of the intersite spin correlations, 
so that (Si • S 2 ) and (Si • S 2 )o almost overlap. 

Although our results presented here are only a starting point and a systematic 
analysis of the properties of the ground states and spin-orbital correlations in realistic 
two-band models with e g and t 2g orbitals of higher dimension is required, our study shows 
that the description of transition metal oxides with partly filled (almost) degenerate 
orbitals has to involve correct symmetry of the orbital degrees of freedom, and thus 
has to go beyond the simplified Hubbard model with degenerate equivalent orbitals. 
This observation is also supported by our recent study of stripe phases in the nickelates 
|46j . Indeed, the diagonal stripe structures with filling of nearly one hole per atom, as 
observed experimentally, are the ground state of the model with the physically relevant 
hopping elements between e g orbitals, while instead the most stable stripes are half- 
filled within the doubly degenerate Hubbard model. We therefore conclude that the 
ground state properties strongly depend on the orbital degrees of freedom active in a 
given compound. 
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Appendix A. Construction of the basis for the two-site model 

We construct a basis of the Hilbert space, starting with the S z = 1 subspace. There are 



two states with T"° = ±1, 

1%) = 4a4al°>> I**t> = ct i €T 4al°>' (A - X) 

and four with T 1 * = 0, 

l^i?) = 775 ( c i^T c 2CT c ict c 2§t)|0)' (A-2) 

l^2|) = 7/2 ^ c i€T c iCT c 2^t c 2Ct)I^)' ( A -3) 
In the S z = subspace there are eight states with — ±1, 

I^) = 72( ct ia4a ±ct ia c 2a)l°)' ( A - 4 ) 

l $ 2 ) = 77 3 (4 a c t ia ±c 2a c 2a )|0), (A.5) 

l $ 3 ) = j=Ma c m ± 4 5T 4 a )l°>> ( A - 6 ) 

|*J> = 775(4 ?T 4a ± 4 a 4 a )l°)' ( A - 7 ) 

and eight with = 0, 

1^5") = \ (4^t4ci + 4^i4ct (4^t4cj. ~^ 4^i4ct)) l^)> ( A -^) 

l^) = I (4^t4ci + 4^i4ct (4a4?t + 4ct4?j.)J |o)> ( A -9) 

i$ 7 ) = 2 ^4^t4ci ~~ 4^i4ct ^ (4ft 4ci — 4fj.4ct)) (a.io) 

1^8 ) = 2 (4ft4cj. — 4fj.4ct ( c ici c 2ft — 4ct4§i)j ( A -n) 



The states |$+) and |$4 ), together with |$ 5 ) and |$g ) belong to the triplet subspace. 

Appendix B. Eigenvalues of the Hamiltonian (JTJ) with t 2g orbitals 

Below we give a complete list of the eigenvalues of Hamiltonian ((TJ for ti g orbitals in 
the absence of crystal field splitting (at E = 0) with the degeneracy of the states given 
in parenthesis: 

— S = 1 subspace 

\\U - 3J H ± y/(U - 3J H ) 2 + 16t 2 } (3), U-3J H (3), 0(9); (B.l) 
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and 



S = subspace 

\\U -J H ± y/(U - J H ) 2 + 16t 2 } (2), U-J H (2), (B.2) 
±{[7 + J H ± v/(C/ + J ff ) 2 + 16t 2 }, £/ + Jjr, 0. (B.3) 



Appendix C. Eigenvalues of the Hamiltonian (pQ) with e 3 orbitals 

Finally we present the eigenvalues of the Hamiltonian given by equation for e g 
orbitals in the absence of crystal field splitting (at E = 0). The degeneracy of the 
states is given in parenthesis: 

— S = 1 subspace 

\\U - 3J H ± y/(U - 3J H y + 4t 2 } (6), 0(6); (C.l) 

— S = subspace 

\[U- J H ± y/{U-J H ) 2 + 4t*} (2), U±J H , 0, (C.2) 

and, 



2U f I rJu\ 2 TTrtV 2 sir -at 



with 



* -t{ 1 + V 1+ »(t) + -ffl}. < C5 » 



1 - (2^) 2 + 2(|f) 2 
cos (a) = - 3-. (^-6) 



+ 3(f) 2 + 12(A) 2 | 2 
The latter three eigenvalues follow from the submatrix in equation (|12jl. 
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